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ABSTRACT 



In the standard scenario of the Newtonian gravity, a late- type galaxy (i.e., a 
spiral galaxy) is well described by a disk and a bulge embedded in a halo mainly 
composed by dark matter. In Brownstein-Moffat gravity, there is a claim that 
late-type galaxy systems would not need to have halos, avoiding as a result the 
dark matter problem, i.e., a modified gravity (non- Newtonian) would account for 
the galactic structure with no need of dark matter. In the present paper, we probe 
this claim via numerical simulations. Instead of using a "static galaxy," where 
the centrifugal equilibrium is usually adopted, we probe the Brownstein-Moffat 
gravity dynamically via numerical iV-body simulations. 

Subject headings: galaxies: spiral - gravitation - methods: n-body simulations 
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Introduction 



Spiral galaxies are very common in the universe. When studied with spectroscopic 
methods, these beautiful objects usually reveal to rotate faster than it would do if only 
their visible matter were taken into account. This is the missing mass problem. Then, one 
generally believes that there is a dark matter component surrounding their disks, making 
them to rotate in that way. In other words, the canonical model says that spiral galaxy 
rotation curves inferred from observations could be accounted for the dark matter halos that 
embed the galaxies (see, e.g., 



Sofue fc Rubin 



2001 



van Albada et al 



19851 ). Therefore, this 



additional distribution of mass causes an additional acceleration and the rotation curves 
can be accounted for. 

This picture to the understanding of the structure of late-type morphologies finds 
support in the present canonical physics, that considers the general relativity the most 
reliable theory of gravitation, and, as a consequence, due to the correspondence principle, 
the Newtonian gravitation at the weak field limit. 

However, in order to avoid the dark matter problem, many authors consider 
modifications in the Newtonian gravity at large scales, in other words, they consider that, 
in galactic or even cosmological scales, the law of gravity is different from the Newtonian 
one. In this way, contrary to the dark matter paradigm, these authors argue that there is 



only a disk, but due to another kind of law of gravity at large scales, it could we 



1 rotate 



in accord with the observed rotat 


ion curves (see, e.g., 


Vlilerom & Sanders 


2005: 


Milgrom 


2008: 


Moffat & Toth 


2008. 


2009: 


Brownstein & Moffat 


2006a 


, and references therein). 









It is worth mentioning that 



Brownstein &: Moffat 



( I2006bl ) also argue that it is possible 



to explain the X-ray cluster masses without dark matter using this very theory. 
In particular, the present work is concerned with the study made by 



Brownstein fc Moffat 
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( I2006al hereafter BM), who present an alternative approach to explain the rotation curv es 



of spir al galaxies, using an alternative theory of gravitation studied by 



Moffat k Toth 



( 12008 



, hereafter Moffatian gravity). Our main aim is to probe this theory using numerical 



simulations. 

The Moffatian gravity is a covariant theory of gravity based on the coupling of the 
Einsteinian gravity with a massive skew symmetric field F a p 7 . Here, we do not present and 
discuss this alternative theory of gravity in detail, since we are only concerned with its weak 
field limit. 

BM claim that in the weak field limit the gravitational potential of a point source of 
mass M reads 
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where G is the Newtonian constant of gravity, r is the distance to the source, M is 
a coupling parameter and tq is a characteristic length. From the BM potential, the 
gravitational acceleration by a particle of mass M reads 
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BM state that the above equations hold in any environment. To apply these equations 
to model their rotation curves, they consider the following procedure: they integrate 
these equations for a spherically symmetric system. As a result, they obtain the following 
equations: 
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(3) 
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and 



a(r) 



GM(r) 



f \ I f 
1 - exp 1 + — 



(4) 



Now, the M outside the braces becomes M(r), but the M inside the braces does not 
change. Therefore, for given values of ro, Mo and M (the total mass of the system) the 
expression inside the braces depends only on r. Assuming centrifugal equilibrium and a 
given mass distribution, i.e. M(r), BM model their rotation curves. 

Also, note that Equation (j4]) obeys the superposition principle, which is an essential 
property if one intends to make use of iV-body simulations. 

We note from the above equations that the Moffatian gravity can be adjusted to a large 
range of lengths and masses, by means of the parameters r and M , that depend on the 
system under investigation. Hence, they could, in principle, explain many physical effects 
that the canonical theory attributes to the dark matter: galaxy rotation curves, galactic 
cluster masses and so on, because the BM's parameters depend on the total mass of the 
system under study. 

BM apply their modified gravity to explain the rotation curves of a sample of observed 



spiral galaxies a nd yield the following set o: 



r = 13.96 kpc ( IBrownstein fc Moffat 



parameters, namely, M = 96 x 1O 1O M and 



2006a ) . As usual in this kind of study, they consider 



centrifugal equilibrium, which means in the end that their galaxy model is "static." 

It is well known, from the iV-body simulation studies, that disks built with exponential- 
Spitzer laws (hereafter SdMH disks), when submitted to the Newtonian potentials, reveal 



to have secular centrifugal equilibrium, 



simulated time (jSpringel &: White 



1999 



maintaining their density profiles 



Springel. Di Matteo fc Hernquist 



or gi 



gayears of 



2005|)- However, 



under the BM's picture, we do not know how long the BM potential can maintain galactic 
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disks in dynamical equilibrium, and in fact if it really can. This is the main aim of the 
present paper. 



In previous works (IBrandao fc de Araujo ll2010al Jblid). we developed an iV-body method 
to study alternative theories of gravity at galactic scales. In particular, we have studied an 
Yukawian gravitational potential and have shown that this potential is viable only if the 
Yukawian parameter is such that the Yukawian potential is nearly Newtonian. 

In the present work, we also use iV-body simulations to construct and evolve spiral 
galaxies now submitted to the potential and acceleration given by Equations (131) and 
PJ, respectively. We adopted, in these simulations, the same set of fitted parameters 
(M = 96 x 1O 1O M and r = 13.96 kpc) as by BM. 

Our aim is to probe if the BM model is dynamically consistent, as it should be in order 
to be considered a realistic model. In this way, we are able to test the BM's claim and 
verify the reliability of thi s kind of law of gravity, usi ng the concepts and techniques of the 



galactic dynamics studies (IBinney fc Tremaine 



2008|). 



In Section [2J we present the iV-body code used in this work and the numerical 
techniques used to model spiral galaxies with the BM potential; in Section the results of 
the numerical simulation are present; and finally, in Section HJ we discuss our main findings 
and conclusions. 



2. Numerical Techniques 



2.1. The Modified Code 



We perform ed the simulations with a modified and tested version of the Gadget- 2 code 



(jSpringel 



20051 ). As discussed in our previous papers ( IBrandao fe de Araujo 



2010a 



d), 



-7- 



this choice is based on the collisionless nature of galaxies modeled by particles, and the 
method used by Gadget-2 to build the recursive tree in order to compute forces and 
potentials. 

We have modified the Gadget-2 and replaced the Newtonian potential and acceleration 
by the expressions given by Equations ([3]) and (j3J just in the respective code's instructions. 
We have extensively tested the code efficiency to calculate potentials via the tree method, 
using a typical initial snapshot composed by a disk embedded in a dark matter halo. 

It is worth noting that we included a halo only to test the code accuracy, since in the 
probes of the ability of the BM gravity to generate spiral galaxies with flat rotation curves, 
without dark matter, the halo is not included in the simulations. 

The Newtonian and the Moffatian codes were set up with the tolerance parameter 
9 = 0.8, maximizing the tree code's performance and softening length /^h — 0.33 kpc 
for the dark halo and = 0.15 kpc for the baryonic disk. This set of softening lengths 
yields conservation of the total energy better than other sets. The procedure to choose 
the values for 9, Uh and U ad opted here is analog to that used in our previous papers 



( iBrandao fc de Araujo 



2010al . and references therein). 



2.2. Initial Conditions 

We have studied three models: model I a typical exponential-Spitzer disk with a dark 
matter halo (Newtonian galaxy); model II a second one, derived from the first one, where 
the system is built and the halo is removed before running the simulation with a Moffatian 
potential (pseudo- Moffatian galaxy); and model III a pure BM's disk model (pure Moffatian 
galaxy). Below, we give the reasons why such models are chosen. 



Model I follows the prescriptions described by 



Springel. Di Matteo fe Hernquist 



(12005 
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and references therein), but with the difference that we are now including neither bulges 
nor gas particles, because we do not need a detailed Newtonian model, since our main aim 
is to investigate the dynamical properties of disk galaxies under alternative potentials. The 



reason for this choice has to do with the ir dynamical stability and computationa 



to model them, as explained in detail in 



Springel. Di Matteo &: Hernquist 



simplicity 



teOQfih . Besides, 



it is worth stressing that SdMH disks reproduce all observational features of late-type 
systems. 

I II 1 

In one of our previous papers (IBrandao fc de Araujo Il2010d ). we described the 

principal recipes to model such galaxies. Besides the parameters discussed in the previous 
section, the models are composed by A^aio = 30, 000 particles for the dark matter halo 
and iV disk = 30, 000 for the disk. All these figures yield a conservation of energy that turns 
our galaxy model reliable, as the reader will find later on in this paper. Model I is then a 
Newtonian model for a spiral galaxy, which is consistent with observations, that we use to 
compare with the Moffatian model for spiral galaxies. 

Model II is the same as model I but without dark halo particles, i.e., it is just the 
Model I from it we have extracted the halo particles. It is worth stressing that the halo 
potential is computed to make the rotation curve of the disk of model II as reliable as 
possible-like the observed ones. 

After extracting the halo, the disk should be stable under the BM's potential. This 
procedure is made to test the BM's hypothesis, which claims that disks can exist under 
Moffat's potential without halos, explaining rotation curves without exotic non-baryonic 
dark matter. In this way, we are able to test this hypothesis over a secular evolution of a 
disk. 

Model III is similar to model II, but with the following differences: (1) we do not 
compute the halo's potential and (2) we replace the Newtonian potential originated 
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from the disk (see, e.g. iBinney fc Tremaine 1 120081 . Equation (2.170)) by the BM's 
potential originated from all the positions of disk particles, using the BM's para meters 



M = 96 x 1O 1O M and r = 13.96 kpc, as fitted by 



Brownstein &: Moffat 



fl2006af ) 



In other words, we first build the positions of the particles and then compute their 
pote ntials using Equation (15]) over the mesh, where the velocities will be calculated (see, 



e.g. 



Brandao de Araujo 



2010d ). We use these results to integrate the Jean's equations, 



and then calculate the velocity dispersions. With this procedure, the initial disk is put in 
equilibrium with the BM's potential during the first snapshots of simulated time. We have 
set iVdisk = 30, 000 particles to this disk (hereafter Moffatian disk or BM's disk for short) in 
order to maintain the same resolution set to the Newtonian model. Note that the disk thus 
modeled is not necessarily consistent with observations, as the reader will see later on. 

The following set of parameters is chosen to the SdMH disks (model I)): the total mass 
Af t = v% 00 /(l0GH ) = 0.98 x 1O 12 M , where v 20 o = 160kms _1 is the virial velocity, G is 
the gravitational constant, and Hq = 100 km s -1 Mpc" 1 , is the Hubble constant; the total 
mass of the disk = m^M t , where m& = 0.041 is a dimensionless fraction of the total 

mass; the disk scale length h = 2.74 kpc; the disk vertical scale height z ~ 0.2/t; and the 
spin parameter A = 0.033. We emphasize that this model has mass comparable to the Milk 
Way. 

For model II we also use the same set of parameters used in model I. Recall, however, 
that after building the disk-halo system, we then extract the halo. 

In the BM's disks, we maintain the total disk mass and the exponential-Spitzer 
profiles. We built the disk particles from the density profiles, and used the BM's parameters 
(Af = 96 x 1O 1O M and r = 13.96 kpc). 



In the next section, we present the results of our simulations. 
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Simulations and Results 



3.1. The Newtonian Model 



Although the Newtonian model can be considered as a can onical one by the 



current theory of galactic a stronomy (IBinney &: Merrifield 



(IBinney &: Tremaine 



19981 ) and galactic dynamics 



2008), it is important to bear in mind that here the Newtonian model 



is not the main result of our work. We have run the Newtonian simulations to compare it 



with models II and III. As it is we 



(jSpringel. Di Matteo fc Hernquist 



1 known from simulations of canoni cal late-type systems 



2005 



Brandao fc de Araujo Il2010d ). many observational 



features, as for example, the exponential radial density profile, the rotation curves, etc can 
be accounted for Newtonian models. 

In Figures (DO)-®, we display the principal results of the Newtonian model we 
simulated. Figure [1] displays a picture with four boxes: at the top left, phase space points 
rxwof the initial snapshot, where r is the total distance from origin (0,0,0) in kpc and v 
is the modulus of the velocity in kms' 1 ; at the top right, the relative energy conservation 
is AE/Eq, where AE = E(t) — E(0) and E = E(0), where E(t) denotes the total energy 
of the system at time t; at bottom left, the phase space points at final time t — 1 Gyr; 
and at the bottom right, the rotation curve is R x v r , where R is the radial cylindrical 
coordinate and v r is the rotation velocity. Note the similarity between the initial and the 1 
Gyr rotation curves; one concludes that early-type galaxy structure can be accounted for 
the Newtonian model quite well, as is well known. Later on we will see that the comparison 
between the density profiles at these same different times corroborates such a conclusion. 

Figures [2] and [3] show the particle's positions in the ^-projection at respective times 
indicated in the boxes. We note from these figures that the simulation follows a typical 
behavior with the formation of spiral arms and central bar. These features are expected 
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Brandao fc de Araujo 



fl2010d . 



from simulations of late-type morphologies, as explained by 
and references therein). In the first 0.33 Gyr, the disk has spiral arms and the system 
evolves to a central bar with two spiral arms afterward. In other words, this simulation 
is reliable and represents a standard typical late-type system without external influences. 
Energy conservation is better than 0.1%. 



3.2. Moffatian Models 

We show above that the Newtonian disks maintain their physical properties for many 
crossing times. Now, when a spiral galaxy is submitted to the BM potential, an anomalous 
behavior occurs to the disk's structure, as we will see in this section. 

In Figure HI we display the same as in Figure [TJ In the top right panel of Figure HI we 
display the evolution of the energy violation. Note that after 1 Gyr (end of our simulation), 
the violation is less than 1%, showing the quality of our simulation. 

From this simulation, we also see that the particles in the phase space are more 
scattered as compared to the Newtonian simulation. This scattering is interpreted as a loss 
of information from the initial configuration. Note, at the bottom right panel, that the 
disk rotates slower than in its initial time, even though it still resembles a typical rotation 
curve. At first glance, it seems, by only examining Figure HJ that the Moffatian potential 
can really explain BM's disks, and no dark matter is needed to maintain secularly stable 
exponential disks. 

We display snapshots of the disk's particles in Figures and El Examining these 
snapshots, we see that at 0.08, 0.16, 0.24, and 0.32 Gyr times the system presents atypical 
substructures, in particular in its central region (R < 10.0 kpc). As the disk evolves, an 
annular substructure appears, and the central region becomes almost empty of particles, a 
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void. Subsequently, the big ring collapses and fragments into three spiral unbound arms. 
At this stage, there is no core. 

At t = 0.66 Gyr, the arms merge themselves and become a bar-like structure surrounded 
by two weak spiral arms. This bar-like structure is small when compared to the Newtonian 
case. This anomalous behavior is due to the fact that the particles are not in perfect initial 
equilibrium with the Moffatian potential. 

However, it is important to bear in mind that we are testing the BM claim, which 
considers that the Moffatian potential can explain many observed rotation curves of galactic 
disks without dark matter halos. We recall that our simulated initial disk reproduces 
the observational features of real galaxies: exponential-Spitzer disk, velocity dispersions, 
rotation curves and baryonic mass, and has secular equilibrium in the Newtonian case. 
From this point of view, we would expect that, if BM hypotheses were reliable, our disk 
would be stable all the time in the simulation. We obtain that an equilibrium configuration 
is reached only after t > 0.66 Gyr. After this time, the system is similar, in some dynamical 
aspects, to the Newtonian case, but we will see below that this similarity is only apparent. 

The results of model II's simulations bring some interesting questions. One could claim 
that the simulations of model II need a fine tuning of the parameters M and r to maintain 
the entire disk in equilibrium. These parameters could, therefore, be changed to avoid 
atypical rings and central voids during the first 0.33 Gyr. To study these sub-morphology 
features, we have made some tests changing the parameters Mq and tq, but we have not 
found any set of parameters that left the simulations just like the Newtonian one, which, it 
is worth stressing again, is consistent with observations. 

It was just the results of model II that lead us to make a pure Moffatian disk; our 
model III whose simulations were made with the same parameters used in the models I and 
II. The results of model III can be seen in Figures [7H9J 
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In the top right panel of Figure [TJ one finds the phase-space points. These points are 
distributed somewhat differently as compared to previous models, because the rotation 
curve has now a different shape. Also, after 1 Gyr of simulated time, we note that the 
scattering in the phase space plane is more prominent than in the model II and the initial 
configuration is completely lost. 

This effect is neither due to simulation inaccuracies nor to the resolution of the 
simulation, as can be seen that the energy violation (see the top right panel of Figure [7]) is 
small, namely, logAE/E < —2.0. 

Our results show in the end that the exponential profile is not consistent with the 
Moffatian potential. As occurred with model II, the rotation curve of the model III is below 
the initial one. For R < 10 kpc, the rotation curve is linear and associated with the central 
core movement, similar to a rigid body rotation. For R > 10 kpc, the rotation curve is 
similar to a typical rotation curve of late-type systems, but, in analogy with model II, this 
similarity is only apparent, as we will see below. 

Figures [8] and [9] display the particle's positions at (x, y) plane, where the simulated 
time is indicated in the respective boxes. Note that, at initial times, the same atypical 
central void appears with its surrounding ring, that collapses and fragments in two small 
clumps, that evolves separately eventually merging at t ~ 0.9 Gyr. 

This simulation shows that for a typical initial equilibrium configuration, BM's disk 
systems evolve in such way that we cannot produce models in agreement with the observed 
values, since atypical configurations appear. 

The mixing in the phase space, for example, shows us that our BM's models do not 
recover the observational characteristics and, therefore, this model cannot be considered as 
a good representation of disk galaxies. 
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But, one could argue that the snapshot at time t — 1 Gyr resembles a typical 
late-type system, and this model could be considered reliable, since after t > 0.9 Gyr this 
configuration seems to be stable. Later on we consider this issue again by studying the 
density profile of the disk. 

We now compare the rotation curve that comes out from our simulations with the one 
that comes from the centrifugal equilibrium, i.e. v 2 (r) = ra(r), where a(r) is given by 
Equation (j4j). This rotation curve reads 

<r) = /^fS {l + /f[l- ex P (-r/r )(l + r/r )] J \ (5) 

where M(r) is obtained from our simulation. Note that the above equation is just the one 
used by BM. 

In Figure [101 we show the result of this comparison. Note the similarity between the 
rotation curves, even when considered different simulated times. The differences come from 
particle noise, very common in iV-body samplings after evolving the system. This result 
shows the consistency between our simulations and the calculation done by BM. 

An important issue is to examine the radial density profile of the simulated disks. A 
reliable model of a spiral galaxy must have at the end of the simulation the same initial 
profile. 

In Figure [TTJ, the radial density profile of the simulated disks is calculated at t — 
and t = 1 Gyr. This figure presents the logarithmic of particle counts per unit of area, 
i.e., log(p), in the (x,y) plane, that was divided in many concentric rings, in cylindrical 
coordinate system. We have also estimated the initial exponential profile analytically, which 
is indicated by dotted lines. 
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Frame I indicates the Newtonian initial profile, estimated for all the particles and 
analytically both profiles coincide exactly. Frame II shows the final Newtonian profile; 
Frame III shows the final profile of the model II; and Frame IV shows the final BM profile 
(model III). Each frame has its initial exponential profile for comparison. Note that the 
Newtonian disk maintains its initial profile, in spite of particle noise, very common in 
iV-body samplings after evolving the system. This figure shows that bars and spiral arms 
left the exponential profile almost unchanged. 

For the Frame III, however, we conclude that the simulations of the model II have 
minimized the central density and maximized the density at R > 10 kpc. In the phase 
space, this can be noted by the migration of particles to more distant regions, where their 
orbits are more stable, increasing the density there. The same happens with BM disks, as 
we can observe in Frame IV. 

It is worth mentioning again that our BM disk has the same initial profile of the 
Newtonian model, and Frame I represents all initial profiles, due to the fact that all initial 
disks are modeled by the same exponential-Spitzer law. We note that, in spite of the final 
snapshot configurations of the models II and III mimic disks, their density profiles are 
very different from a canonical exponential-Spitzer law, from which they were built. As a 
conclusion, the Moffatian potential cannot maintain exponential-Spitzer disks in dynamical 
equilibrium and therefore cannot be consistent with observations. 

In Figure [121 we present, for model III, the radial density profile for different simulated 
times. This figure clearly shows that, even starting from a disk in equilibrium, the Moffatian 
gravity does not maintain the exponential-Spitzer disk with the same initial density profile. 
As time goes on, the profile turns gradually and is definitely different from the exponential 
profile. This shows that simulating models are a more robust tool than simply adjust 
profiles in static models, because the evolution of the system cannot be followed in this last 
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approach. 

In conclusion, an alternative gravitational law must explain not only the rotation 
curves of spiral galaxies, but also their density profiles. In this way, simulations can be a 
powerful tool to deal with such an issue. 

3.2.1. An additional test 

Note that depending on the values of and Mo, the Newtonian gravitation would be 
naturally recovered. Therefore, as an additional test we examine such an issue, which could 
help one to see how the structure of a disk galaxy would be modified by an alternative 
gravity law. 

We start from a Moffatian disk, just the one present, for example, in the first snapshot 
of Figure El whose initial rotation curve is shown in Figure [7] Then, we set r = lOOOkpc 
and Mo = 1 x 10 10 M and follow how the disk evolves. 

The results of this simulation are shown in Figure [13] and HH For comparison, we 
also performed simulations using the Newtonian Gadget-2. Both the calculations produce 
identical results. Therefore, the Newtonian calculation is recovered. 

In particular, note in Figure US] the rotation curve after 1 Gyr of simulated time; for 
large values of R, the velocity is a decreasing function of R, as expected for a Newtonian 
disk. 

It is also interesting to compare the Moffatian disk snapshots and Newtonian ones. 
Comparing Figures [8] and HU one clearly sees how different these configurations are. 
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4. Conclusions 

Studies of alternative theories of gravity are usually made statically, instead of modeling 
"live" systems. For example, in the study of spiral galaxies centrifugal equilibrium is usually 
considered. The problem with this approach is that the secular evolution of these systems 
cannot be followed. 



In the present paper, in particular, we follow the techniques used by lBrandao fc de Araujo 



(12010aU b 



cj) and probe the BM's model, verifying if the Moffatian gravity is dynamically 
consistent to account for spiral galaxies. A dynamical study is a robust approach, since 
it considers the distribution functions and Jean's equations that can maintain A^-body 
systems in equilibrium. 

We performed the simulations with a modified and tested version of the Gadget-2 
code. We replaced the Newtonian potential and acceleration by the expressions given 
by Equations (ED and (j4j) just in the respective code's instructions. We have extensively 
tested the code efficiency to calculate potentials via the tree method, using a typical initial 
snapshot composed by a disk embedded in a dark matter halo. 

From our simulations, it follows that the Moffatian potential cannot generate 
exponential disks in equilibrium, although the rotation curves are nearly flat throughout 
the whole simulated time. The final (stable) equilibrium configuration is very different from 
the exponential profile, as we have explained above. So, if we consider that exponential 
profiles are reliable, as it seems to be the case, since they are consistent with observations, 
it is hard to believe that BM model is better than the Newtonian model (baryonic disk 
embedded in a dark matter h alo) to account for the spiral galaxy structure. Recall that 



Brownstein fc Moffat 



(l2006al ) claim that M and r fit very well many (observed) rotation 



curves, but they have only analyzed static models (centrifugal equilibrium configurations). 
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Finally, it is worth mentioning that it would be very interesting to simulate the 
BM potential for standard elliptical systems and to probe its influence over their global 
properties, like we have done in this work. A forthcoming paper will consider such an issue. 

We are grateful to the Brazilian agencies CAPES, CNPq and Fapesp for support. We 
also thank Nick Gnedin for providing us a visualization software. Last but not least, we 
thank the referee for useful suggestions and criticisms, which greatly improved the final 
version of our paper. 
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Fig. 1. — Top left: phase space for initial snapshot data. Top right: energy conservation of 
the simulation. Bottom left: phase space for final snapshot data at 1 Gyr. Bottom right: 
rotation curves for initial (solid line) and final (dashed line) snapshots. 
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Fig. 2. — Newtonian disk at ^-projection at 0, 0.33, 0.66, and 1 Gyr of simulated time 
(indicated in the respective boxes). 
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Fig. 3. — First 300 Myr of simulated time to the Newtonian disk at ^-projection. Time is 
indicated in the respective boxes. 
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Fig. 4. — Top left: phase space for initial snapshot data. Top right: energy conservation of 
the simulation. Bottom left: phase space for final snapshot data at 1 Gyr. Bottom right: 
rotation curves for initial (solid line) and final (dashed line) snapshots. 
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Fig. 5. — Pseudo-Newtonian disk at z-projection at 0, 0.33, 0.66 and 1 Gyr of simulated 
time (indicated in the respective boxes). 
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Fig. 6. — First 300 Myr of simulated time to the pseudo-Newtonian disk at ^-projection. 
Time is indicated in the respective boxes. 
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Fig. 7. — Top left: phase space for initial snapshot data. Top right: energy conservation of 
the simulation. Bottom left: phase space for final snapshot data at 1 Gyr. Bottom right: 
rotation curves for initial (solid line) and final (dashed line) snapshots. 
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Fig. 8. — BM disk at ^=projection at 0, 0.33, 0.66, and 1 Gyr of simulated time (indicated 
in the respective boxes). 
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Fig. 9. — First 300 Myr of simulated time to the MB disk at z-projection. Time is indicated 
in the respective boxes. 
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Fig. 10. — Moffatian rotation curves (model III). Solid lines stand for the rotation curve 
given by Equation (jSJ) and dashed lines for the rotation curves obtained from the snapshots. 
Time is indicated in the respective boxes. 
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Fig. 11. — Radial density profiles of the simulated models. Horizontal axis: radial distance 
R from the center of gravity. Vertical axis shows the log of the particle counts, log(p), where 
p is the number of particles per area unit. Frame I shows the initial values of the density 
profile, represented by solid lines. Frame II: the final profile of the Newtonian model. Frame 
III: the profile of the pseudo-Moffatian disk, and Frame IV: the profile of the Moffatian disk. 
All frames display the exponential profile in dotted lines that match with the initial counts. 
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Fig. 12. — Evolution of the radial density profiles of the model III. Time is indicated in the 
respective boxes. 
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Fig. 13. — Model III with the same initial snapshot used in Figure obtained with the 
Gadget-2 modified code in the Newtonian limit, i.e., with r = 1000 kpc and M = 1 x 
1O 1O M . Top left: phase space for initial snapshot data. Top right: energy conservation of 
the simulation. Bottom left: phase space for final snapshot data at 1 Gyr. Bottom right: 
rotation curves for initial (solid line) and final (dashed line) snapshots. 
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Fig. 14. — BM disk at z-projection at 0, 0.33, 0.66, and 1 Gyr of simulated time (indicated 
in the respective boxes) obtained with the Gadget-2 modified code in the Newtonian limit, 
i.e., with r = 1000 kpc and M = 1 x 1O 1O M . 



